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We present the first spin alignment measurements for the K*^{9Q2) and (j!>(1020) vector mesons 
produced at mid-rapidity with transverse momenta up to 5GeV/c at ^/snn = 200 GeV at RHIC. 
The diagonal spin density matrix elements with respect to the reaction plane in Au+Au collisions 
are poo = 0.32 ± 0.04 (stat) ± 0.09 (syst) for the K*° (0.8 < pr < 5.0 GeV/c) and poo = 0.34 ± 0.02 
(stat) ± 0.03 (syst) for the (f) (0.4 < pr < 5.0 GeV/c), and are constant with transverse momentum 
and collision centrality. The data are consistent with the unpolarized expectation of 1/3 and thus 
no evidence is found for the transfer of the orbital angular momentum of the colliding system to the 
vector meson spins. Spin alignments for A"*" and (f) in Au+Au collisions were also measured with 
respect to the particle's production plane. The (j> result, poo ~ 0.41 ± 0.02 (stat) ± 0.04 (syst), 
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is consistent with that in p+p colhsions, poo = 0.39 ± 0.03 (stat) ± 0.06 (syst), also measured in 
this work. The measurements thus constrain the possible size of polarization phenomena in the 
production dynamics of vector mesons. 

PACS numbers: 24.70.+S, 25.75.Nq 



Measurements of the polarization of the particles pro- 
duced in relativistic heavy-ion collisions may provide new 
insights into the initial conditions and evolution of the 
nuclear system during the collision [J, [3, 9 ■ particu- 
lar, by studying the polarization magnitudes with respect 
to various kinematic planes one could attempt to discern 
the point in the evolution of the system at which the po- 
larization arises and, hence, the dominant mechanisms 
involved. The planes that are relevant to this paper are 
the reaction plane, which is defined by the beam momen- 
tum and the nuclear impact parameter, and the particle's 
production plane, which is defined by the beam momen- 
tum and the momentum of the produced particle. 

In non-central relativistic heavy-ion collisions, trans- 
verse gradients of the total longitudinal momentum of 
the participant matter result in substantial local orbital 
angular momentum of the created partons [l| . Due to the 
spin-orbit coupling of QCD, this orbital motion may re- 
sult in a net polarization of the produced particles along 
the direction of the initial angular momentum, that is, 
perpendicular to the reaction plane, yielding a global 
hadronic polarization in the final state |1|,|J,|5|]. The mag- 
nitude and the transverse momentum (px) dependence of 
the global polarization are therefore expected to be sensi- 
tive to different hadronization scenarios Q . In particular, 
the proposed quark recombination model for hadroniza- 
tion of bulk partonic matter created at RHIC Q , which 
reproduces measurements in the intermediate pr region 
{2 < pt < 5 GcV/c) quite well @, may be an effective 
dynamical mechanism for transferring polarization from 
quarks to vector mesons or hyperons. 

One can study the polarization of final state hadrons 
with respect to the particle's production plane as well. 
Non-zero polarizations transverse to this plane are ex- 
pected to be sensitive to particle formation dynamics 
and topossible intrinsic quark transverse spin distribu- 
tions [9(1 . Because large production plane polarizations 
have been observed for hyperons in unpolarized p-l-p and 
p+A interactions and for vector mesons in K+p, 
n-|-C, and e+e" interactions [O, [13, [Si 13 1 ^he disap- 
pearance of these effects in relativistic heavy-ion colli- 
sions might indicate that the system is isotropic to the 
extent that, locally, there is no longer a preferred direc- 
tion [l^ . It has also been suggested that vector me- 
son spin alignment with respect to the production plane 
is closely related to the single-spin left-right asymmetries 
in transversely polarized p+p collisions [ItI . [TsI . [l9l . [20| . 

We have recently measured the global (reaction plane) 
polarizations of A and A hyperons produced at mid- 
rapidity in Au-|-Au collisions at center-of-mass energies 
y/J^ = 62.4 GeV and 200 GeV [IJ. The results, 
I -Pa A I — 0.02, are well below the predictions of Ref. [l| 



and are in agreement with the refined calculations of 
Ref. 0. 

In this paper, we present the first measurements of the 
spin alignment for the K**^ and (j) vector mesons with re- 
spect to both the reaction plane and the production plane 
in Au+Au collisions at y^sjv77 = 200 GeV. We present 
also the spin alignment of the (f> meson with respect to 
its production plane for p-f p collisions at the same colli- 
sion energy. The data cover 0.8 < pr < 5.0 GeV/c 
and the (j> data cover 0.4 < pT < 5.0 GeV/c. 

Spin alignment is described by a spin density matrix 
p, a 3 X 3 Hcrmitian matrix with unit trace. A devia- 
tion of the diagonal elements Pmrni^n ~ "liO, 1) from 
1/3 signals net spin alignment. Because vector mesons 
decay strongly, the diagonal elements p-i-i and pn are 
degenerate and poo is the only independent observable. 
It can be determined from the angular distribution of the 
decay products [l^ , 

^^^^ = A^o X [(1 - poo) + (3poo - 1) cos2 9*1 (1) 

where Nq is the normalization and 9* is the angle between 
the polarization direction n and the momentum direction 
of a daughter particle in the rest frame of the parent 
vector meson. In the case of a global spin alignment 
measurement, the polarization direction n is along the 
orbital angular momentum of the colliding system. It 
is determined by the reaction plane, requiring only the 
second order term since Eq.(l) is invariant under 9* 
TT + 9* [2^. For the production plane measurement, h 
lies along the normal to the production plane, which is 
determined by the momentum of the vector meson and 
of the colliding beams. 

Vector mesons are expected to originate predominantly 
from primordial production [2^ . [25j , unlike hyperon pro- 
duction which is expected to have large resonance decay 
contributions [2l|, Ull- Another difference between 
the present spin alignment measurement and our recent 
measurement of global hyperon polarization [2]| is that 
contributions to a spin alignment measurement are gener- 
ally additive, whereas contributions along a polarization 
direction may cancel. Last, as far as the reaction plane 
resolution is concerned, the aforementioned method has 
an advantage over the method used in Ref. [2l[, where 
the reaction plane was estimated in forward detectors. 

A total of approximately 2.3x10^ events from Au+Au 
data collection in the year 2004 run and 6.0x10^ events 
from p+p data collection in the year 2001 run have been 
used in these analyses. The events were collected with 
minimum bias triggers [26l . [27j . Charged tracks were re- 
constructed with the STAR Time Projection Chamber 
(TPC) for pseudo-rapidities I77I < 1.0 and all azimuthal 
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angles [28[. Particle identification is achieved by cor- 
relating the ionization energy loss (dE/dx) of charged 
particles in the TPC gas with their measured momenta. 
The measured (dE/dx) is reasonably well described by 
the Bichsel function smeared with a resolution of width 
a [2^. By measuring the (dE/dx), pions and kaons could 
be identified up to a momentum of about 0.6 GcV/c while 
protons could be separated from pions and kaons up to 
a momentum of about 1.1 GeV/c. Tracks within 2(7 of 
the pion/kaon Bichsel curve were selected in the analy- 
ses. The K*'^ and (j> mesons were reconstructed through 
their respective hadronic decay channels, K*^ K ^tt^, 
K-1T+ and K+K'. The K*° and K*° sam- 
ples were combined to enhance the statistics and the term 
i^T*" in the remainder of this paper will refer to the com- 
bined sample. The collision ccntrality was determined 
by the charged hadron multiplicity within \ri\ < 0.5. The 
same analysis techniques have been used in our earlier 
measurements of K*^ and (f> production [2^, [23| . 

Figure [T] illustrates aspects of the data and analysis 
for particular pT bins. The top panels show the invari- 
ant mass distributions for (a) K*^ and (b) (p candidates 
in mid-central (20-60%) Au-f-Au collisions at ^/sj^ = 
200 GeV including all values of cos6'* . In these analyses 
invariant mass distributions were obtained for K*'^ and 
for each cosO* and pt interval. The raw K*'^ and cj) yields 
in each of these distributions were obtained by subtract- 
ing the corresponding combinatorial backgrounds and 
fitting the remaining distributions with a Breit-Wigner 
function plus a polynomial curve to describe the residual 
background. The raw yields were then corrected for de- 
tection efficiency and acceptance determined from Monte 
Carlo GEANT simulations [li,!!^]. The middle panels in 
Fig. [T] show the cos9* distributions, after efficiency and 
acceptance corrections, for the (c) A'*° and (d) (j), respec- 
tively. Eq. ^ was fitted to these distributions to deter- 
mine Poo(pt)- In the analyses, we used charged particle 
tracks with 0.2 < pr < 2.0 GeV/c and pseudo-rapidity 
I77I < 1.0 originating from the primary interaction ver- 
tex to reconstruct the event plane as an estimate of the 
reaction plane [lOl- Tracks associated with a K*^ or a 
(j) candidate are explicitly excluded from the event plane 
calculation. The results for pooiPr) were corrected for 
the finite event plane resolution, which was determined 
by correlating two random sub-events. The correction 
factor on (3poo — 1) is determined to be 1/0.81 in 20-60% 
Au-f Au collisions at ^/snn = 200 GeV [301 ■ The bottom 
panel (e), (f), (g) represents the cos9* distribution for 
K*'^ and mesons with respect to the production plane 
in Au-fAu and p-l-p collisions. In this case, Poo(pt) is 
extracted directly by fitting Eq. ([T|) to the distributions. 
We have checked our analysis procedure by extracting 
Poo for the abundantly produced, but spin-less Kg me- 
son (J^ = 0"). The results are shown in panel (h) and 
are consistent with 1/3 within the statistical uncertain- 
ties, as expected. The x^/nd/ value is unsatisfactory for 
the isTg fit for 0.4 < px < 0.8 GeV/c, which is indica- 
tive of point-to-point systcmatics. It reaches satisfactory 
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FIG. 1: (color online) The invariant mass distribution after 
combinatorial background subtraction for (a) the and (b) 
the <j) meson in mid-central Au-fAu collisions at ^^snn= 200 
GeV including all values of cos.9* . The continuous lines rep- 
resent fits of signal, described with a Breit-Wigner function, 
and residual background, described with the dashed polyno- 
mial curves. Panel (c), (e) and panel (d), (f) represent the 
cos9* distributions for the and </> yields in Au-|-Au colli- 
sions, respectively. Panel (g) is the 4> yield in p-l-p collisions 
while panel (h) shows the control measurement of the spin-less 
A's meson cos9* distribution. The error bars show statistical 
uncertainties. The blue dashed lines in (c)-(h) are fits of Eq. 
(1) to the data. 
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TABLE I: The averaged spin density matrix elements poo with 
respect to the reaction plane in mid-central Au+Au collisions 
at ^/s^ = 200 GeV. 
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FIG. 2: (color online) The spin density matrix elements poo 
with respect to the reaction plane in mid-central Au-|-Au col- 
lisions at y^sjvjv ~ 200 GeV versus pr of the vector meson. 
The sizes of the statistical uncertainties are indicated by er- 
ror bars, and the systematic uncertainties by caps . The A"*" 
data points have been shifted slightly in pr for clarity. The 
dashed horizontal line indicates the unpolarized expectation 
Poo = 1/3. The bands and continuous horizontal lines show 
predictions discussed in the text. 

The measurements of the A'*° and (/> global spin align- 
ment versus px of the vector meson for mid-central 
Au-|-Au collisions at ^/sNN = 200 GeV arc presented 
in Fig. [21 and are summarized in Table |T1 Both statis- 
tical and systematic uncertainties are shown. The dom- 
inant contribution to the systematic uncertainty for the 
(/) {K*°) meson ranges from 0.020 (0.05) to 0.045 (0.10), 
originating from uncertainty in the magnitude and shape 
of the residual background after the subtraction of com- 
binatorial background. This residual arises from the in- 
complete description of combinatorial background via the 
event mixing technique and from distortions to the back- 
ground in the invariant mass distribution near the peak 
caused by photon conversions and other correlated back- 
grounds such as at"* — !■ AT+TT", p° tt+tt", A — *■ pir^ 
and A Nn decays [3l|. In the case of the A'*° 
these backgrounds include A"^ tt+tt^, p° — tt^tt", 
(f> K'^K^, A p7r~ and A Nn decays [3^ . 
Other point-to-point systematic uncertainty associated 
with particle identification for the (AT*") meson were 



estimated to range from 0.007 (0.06) to 0.012 (0.09) by 
tightening the (tt and AT) {dE/dx) cut from 2a to la. 
An additional sizable contribution to the (/) uncertainty 
was estimated to range from 0.007 to 0.012 by varying 
the fitted invariant mass range from 1.00-1.04 GeV/c^ 
to 1.00-1.06 GeV/c^, and to the A'*" uncertainty rang- 
ing from 0.02 to 0.05 by changing its analyzed rapidity 
range from |y| < 1 to \y\ < 0.5. The systematic uncer- 
tainties in the K^* measurements are larger than those 
in the (j) measurement mainly because of the lower signal- 
to-background ratio of ^ 1/1000 compared to ~ 1/25 for 
the (j) meson. The contributions to the systematic uncer- 
tainty caused by elliptic flow effects and the event plane 
resolution are found to be negligible. The Ar*° and (j) 
data are consistent with each other and are consistent 
with 1/3 at all px- 

Hadronization of globally polarized thermal quarks, 
typically having px < IGeV/c, in mid-central Au+Au 
collisions is predicted to cause px dependent deviations 
of Poo from the unpolarized value of 1 /3 [l|, 0, 0, [13 ■ Re- 
combination of polarized thermal quarks and anti-quarks 
is expected to dominate for px < 2GeV/c and leads to 
values of poo < 1/3 as indicated in Fig. [2] for a typi- 
cal range of expected light (strange) quark polarizations 
Pq{s) M- The fragmentation of polarized thermal quarks 
with larger px, however, would lead to values of poo 
> 1/3 for 1 < pT < 3 GeV/c 0, H, which is indi- 
cated as well. In the region of 1 < < 2 GeV/c both 
hadronization mechanisms could occur and their effects 
on Poo may cancel. As observed in Fig. [2] these effects 
are predicted to be smaller than our experiment sensi- 
tivity. However, the large (strange) quark polarization, 
Pq s = —0.3, considered in the recombination scenario 
of Ref. [l| results in worse agreement of poo with our 
(j) data than -0.03 < Pq,s < 0.15 discussed in Ref. [1]. 
Our data are consistent with the unpolarized expectation 
Poo = 1/3. Recent measurement of the A and A global 
polarization also found no significant polarization and an 
upper limit, \Pjyj^\ < 0.02, was obtained [2ll |. 

The centrality dependence of the global spin align- 
ment measurements for K*'^ and (j) vector mesons with 
low and intermediate px is shown in Fig. [3l The or- 
bital angular momentum of the colliding system depends 
strongly on the collision centrality. Global polarization 
is predicted to be vanishingly small in central collisions 
and to increase almost linearly with impact parameter 
in semi-central collisions due to increasing particle angu- 
lar momentum along with effects of spin-orbit coupling 
m QCD Q. The data exhibit no significant spin align- 
ment at any collision centrality and thus can constrain 
the possible size of spin-orbit couplings. 

Figurc|4]and Table Ull present the A'*° and <j> spin align- 
ment measurements with respect to the production plane 
in mid-central Au-|-Au collisions at ^/SNN = 200 GeV to- 
gether with the cf> meson results in p-l-p collisions at the 
same incident energy. As is the case for our measure- 
ments with respect to the reaction plane, the uncertain- 
tics in the measurement with respect to the production 
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FIG. 3: (color online) The dependence of poo with respect 
to the reaction plane on the number of participants at mid- 
rapidity in Au+Au collisions at ^/sNN = 200 GeV. The sizes 
of the statistical uncertainties are indicated by error bars, and 
the systematic uncertainties by caps. The data for pr > 
2 GeV/c and the Jf*" data points have been shifted slightly 
in {Npart} for clarity. The dashed horizontal line indicates the 
unpolarized expectation poo — 1/3. 
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FIG. 4: (color online) The spin density matrix elements poo 
with respect to the production plane in mid-central Au-|-Au 
and p+p collisions at ^/sNN = 200 GeV versus pr of the 
vector meson. The sizes of the statistical uncertainties are 
indicated by error bars, and the systematic uncertainties by 
caps. The Jf*" and the (fi p+p data points have been shifted 
slightly in pr for clarity. The dashed horizontal line indicates 
the unpolarized expectation poo = 1/3. 



plane are smaller for the if) than for the A'*" meson, and 
the statistieal uncertainties are somewhat smaller than 
the systematic uncertainty estimates. The </> point-to- 
point systematic uncertainty estimate includes a dom- 
inant contribution ranging from 0.030 to 0.045 due to 
residual background plus two smaller contributions of 
0.006-0.012 and 0.005-0.010 estimated by varying the 
identification cut on (dE/dx) from 2a to la and 
the fit range of the meson invariant mass from 1.00- 
1.04 GeV/c^ to 1.00-1.06 GeV/c^. For the K*° we esti- 
mate a residual background contribution to the point-to- 
point systematic uncertainty ranging from 0.02 to 0.08 
and about equal contributions ranging from 0.01 to 0.08 
by varying particle identification criteria and analyzed 
rapidity. The Au+Au data for poo a-re consistent with 
1/3 to within 1-2 times the total uncertainties, though 
the central values tend to increase with decreasing pT for 
Pt < 2.0 GeV/c. The p-|-p results are consistent with the 
Au+Au results and with 1/3. No conclusive evidence is 
found for large polarization phenomena in the produc- 
tion dynamics of vector mesons in the covered kinematic 
region with the precision of current measurements. The 
p+p results are in qualitative agreement with the sug- 
gested relation of vector meson spin alignment with re- 
spect to the production plane and the null results ob- 
served for the transverse spin asymmetries in singly po- 
larized p+p collisions at mid-rapidity [l^,!!^]. OPAL and 
DELPHI have previously reported similar null results for 
the spin alignment of the K*'^ and (p mesons produced 
with small fractional momenta {xp < 0.3, Xp = p/pbeam) 
in e+e^ collisions [l^, [l3|, though the production and 
fragmentation processes involved there are different from 
those at RHIC. 



TABLE II: The averaged spin density matrix elements poo 
with respect to the production plane in mid-central Au+Au 
collisions and the </!> result in p+p collisions at ^/sNN = 200 
GeV. 









poo(pT < 2.0 GeV/c) 


0.43 ± 0.04 ± 0.09 


0.42 ± 0.02 ± 0.04 


poo(pT > 2.0 GeV/c) 


0.38 ± 0.04 ± 0.09 


0.38 ± 0.03 ± 0.05 


Poo(PT < 5.0 GeV/c) 


0.42 ± 0.03 ± 0.09 


0.41 ± 0.02 ± 0.04 


poo (p+p) 




0.39 ± 0.03 ± 0.06 



In summary, we have presented the first measurements 
of spin alignment for K*^ and (/) vector mesons at mid- 
rapidity at RHIC. The results for the diagonal spin den- 
sity matrix element poo with respect to the reaction plane 
in Au+Au collisions are found to be constant with pT in 
the measured region, covering up to 5GeV/c, and con- 
stant with centrality. The data are consistent with the 
unpolarized expectation of 1/3 and thus provide no evi- 
dence for global spin alignment despite the large orbital 
angular momentum in non-central Au+Au collisions at 
RHIC. The results with respect to the production plane 
are found to be less than two standard deviations above 
1/3 in Au+Au collisions and arc consistent with the re- 
sults in p+p collisions at the same collisions energy. The 
measurements thus constrain the possible size of polar- 
ization phenomena in the production dynamics of vector 
mesons. Future measurements of polarization with re- 
spect to the jet production plane are complementary to 
the current measurements because they are not sensitive 
to the initial conditions and may probe the system's mean 
free path . 
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